Crossover effects in the bond-diluted Ising model in three dimensions 
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Abstract 

We investigate by Monte Carlo simulations the critical properties of the three-dimensional bond-diluted Ising 
model. The phase diagram is determined by locating the maxima of the magnetic susceptibility and is compared to 
mean-field and effective-medium approximations. The calculation of the size-dependent effective critical exponents 
shows the competition between the different fixed points of the model as a function of the bond dilution. 
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1. Introduction 

The qualitative influence of quenched disorder 
at second-order phase transitions is well under- 
stood since Harris proposed a relevance crite- 
rion [1] based on the knowledge of the specific heat 
critical exponent ct purc of the pure model: when 
a pure is positive, the disordered system will reach 
a new fixed point with new critical exponents 
whereas if a pU re is negative, the same universality 
class will persist. 

As a paradigmatic model, the three-dimensional 
(3D) disordered Ising model characterized by 
ttpurc = 0.109(4) has been extensively studied by: 

• Renormalization group methods in the weak 
quenched dilution regime. The best estimates for 
the critical exponents obtained by this method 
are [2]: 

v = 0.678 ± 0.010, 

r) = 0.030 ± 0.003, 



7 = 1.330 ± 0.017. 

• Monte Carlo simulations of the site-diluted 
case for which the following exponents have been 
found [3]: 

v = 0.6837 ± 0.0053, 
(3 = 0.3546 ± 0.0028, 
7 = 1.342 ± 0.010. 

• Experimental investigations. 

For a review of these different results, see Ref. [4]. 

The general picture which is now widely ac- 
cepted is that, starting from the pure Ising model 
(with exponents v = 0.6304 ± 0.0013, 7 = 
1.2396 ±0.0013, (3 = 0.3258 ± 0.0014, see Ref. [5]), 
the critical temperature decreases and eventually 
vanishes at the percolation threshold below which 
there is no longer any long-range order in the sys- 
tem, due to the absence of percolating clusters. 
In the dilution-temperature plane, the critical line 
T c (p) which delimits ferromagnetic and paramag- 
netic phases has two end-points (namely the pure 
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system and the percolation point) described by 
unstable fixed points and it is commonly believed 
that the quenched disordered system has univer- 
sal properties described by a unique stable fixed 
point. According to this picture, the competition 
between fixed points possibly leads to an effective 
variation of universal critical quantities. 

In this paper, we first determine the phase dia- 
gram of the bond-diluted problem (up to now, only 
the site-diluted case has been studied) and com- 
pare it to a prediction in the single-bond effective- 
medium approximation and then we try to illus- 
trate the competition between fixed points leading 
to crossover regimes in some physical quantities. 



2. Phase diagram 

The bond-diluted Ising model is defined by the 
following Hamiltonian with independent quenched 
random interactions: 



>'// X K ii S "i,"i 
(id) 



{<Ji = ±1). 



(1) 



The coupling strengths are allowed to take two 
different values Kij — K = J/ksT and with 
probabilities p and 1 — p, respectively, 

P(Kij) - pSiKij -K) + (l- p)5(K ij ), (2) 

c = 1 — p being the concentration of missing bonds, 
which play the role of the non-magnetic impu- 
rities. The simulation technique is based on the 
Swendsen-Wang cluster algorithm with periodic 
boundary conditions in the three space directions. 

The phase diagram is obtained numerically from 
the maxima of a diverging quantity (Fig. 1). Here 
we choose the susceptibility, since the stability of 
the disordered fixed point implies that the specific- 
heat exponent is negative in the random system. 
Thus, the error in this quantity is larger than for 
the susceptibility. The percolation threshold is lo- 
cated at p c w 0.2488. 

To get an accurate determination of the max- 
ima of the susceptibility, we used the histogram 
reweighting technique with 2 500 Monte Carlo 
sweeps (MCS) and between 2 500 and 5 000 sam- 
ples of disorder. The number of Monte Carlo 
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Fig. 1. Variation of the average magnetic susceptibility 
Xl versus the coupling strength K = J/kgT for several 
concentrations p and L = 8, 10, 12, 14, 16, 18, 20. For each 
value of p and each size, only one value of K has been 
simulated: the extension of the K values has been obtained 
by the standard histogram reweighting technique. 
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Fig. 2. Variation on a log-log scale of the energy autocor- 
relation time Te versus the size L of the system. The en- 
ergy autocorrelation time increases with the concentration 
of magnetic bonds p and its greatest value obtained for 
L = 96, p = 0.7 is around 9. 

sweeps is justified by the increasing behaviour of 
the energy autocorrelation time and we chose for 
each size at least 250 independent measurements 
of the physical quantities (ATmcs > 250 te)- 

For a second-order phase transition, the auto- 
correlation time is expected to behave as L z at the 
critical point where z is the dynamical critical ex- 
ponent (Fig. 2). For the disordered Ising model, 
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Table 1 

The dynamical critical exponent z as obtained from linear 
fits Of log T£ vs log L. 



V 


1 


0.7 


0.55 


0.4 


z 


0.59 


0.41 


0.38 


0.27 



Table 2 

Evolution of the susceptibility with the number of Monte 
Carlo sweeps per spin for different samples, Xj> an d the 
average value (with 2 500 samples) at L = 96, p = 0.7. 



# MCS 


XI 


X2 


X3 


X4 


XB 


X 


100 


1268 


720 


1141 


939 


833 


1058 


500 


1272 


1520 


1223 


1029 


953 


1210 


1000 


1262 


1544 


1205 


1068 


911 


1219 


1500 


1282 


1433 


1277 


1047 


915 


1227 


2 000 


1332 


1441 


1221 


1073 


917 


1235 


2 500 


1358 


1484 


1234 


1012 


1014 


1234 



we get the values of z shown in Table 1. We see 
that the critical slowing down weakens for the dis- 
ordered model and becomes smaller when the con- 
centration of magnetic bonds p decreases, but it is 
necessary to increase the number of disorder reali- 
sations when p decreases because of the vicinity of 
the percolation threshold. 

In order to check the quality of the averaging 
techniques, we can study the stability of the sus- 
ceptibility for the largest size considered versus the 
number of Monte Carlo sweeps involved in the ther- 
mal average. The results are given in Table 2 for 
different samples as well as for the disorder aver- 
age. With 2 500 MCS, the accuracy of the results 
for a given sample is not perfect, of course, but the 
precision of the average over disorder is quite good 
on the other hand. The disorder average procedure 
has been investigated by computing the suscepti- 
bility Xj f° r different samples, 1 < j < N s , where 
N s is the total number of samples (Fig. 3). We can 
see that the dispersion of the values of x is not very 
large because the fluctuations in the average value 
disappear after a few hundreds realisations. 

The phase diagram obtained from the location of 
the maxima of the susceptibility for the largest lat- 
tice size as a function of the concentration of mag- 
netic bonds is shown in Fig. 4. The simple mean- 
field transition temperature is drawn for compar- 
ison: it gives a linear behaviour as a function of 
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Fig. 3. Distribution of the susceptibility for the different 
disorder realisations of the Ising model with L = 96 and 
a concentration of magnetic bonds p = 0.7. The average 
value over the samples \ is shown by the black line. 
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Fig. 4. Phase diagram of the 3D bond-diluted Ising model 
compared with the mean-field and effective-medium ap- 
proximations. 

the bond concentration p and we can check that 
this approximation holds only in the low-dilution 
regime, p > 0.8. On the other hand, the effective- 
medium approximation [6] gives very good agree- 
ment with the simulated transition line. Treating 
only a single bond in an effective medium leads to 
the following relation for the critical coupling: 



K c {p) = In 



(1 ~Pc)e 



(1-p) 



P-Pc 



(3) 



This relation is exact in the vicinity of both the 
pure system and the percolation threshold. 
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3. Competition between the fixed points 



The main problem encountered in previous stud- 
ies of the disordered Ising model was the question 
of measuring effective or asymptotic exponents. Al- 
though the change of universality class should hap- 
pen theoretically for an arbitrarily low disorder, it 
can be very difficult to measure the new critical ex- 
ponents because the asymptotic behaviour cannot 
always be reached practically. Another difficulty 
comes from the vicinity of the ratios j/v and /3/v 
in the pure and disordered universality classes. In- 
deed, these values for the pure model are [5] : 

7/1/ = f .966(6), /3/u = 0.517(3), v = 0.6304(13), 

and for the disordered Ising model [3] : 

7/1/ = 1.96(3), 0/v = 0.519(8), v = 0.6837(53). 

Thus, from standard finite-size scaling tech- 
niques, the critical exponent v only will allow us 
to discriminate between the two fixed points. This 
exponent can be evaluated from the finite-size 
scaling of the derivative of the magnetisation ver- 
sus the temperature which is expected to behave 
as ~ L x l v . From this power-law behaviour, 

we have extracted the effective size-dependent ex- 
ponent (l/f) c ff which is plotted against 1/L m ; n 
for different bond concentrations p in Fig. 5 where 
L m in is the smallest lattice size used in the fits. 

We clearly see that in the regime of low dilution 
(p close to 1), the system is influenced by the pure 
fixed point. On the other hand, when the bond 
concentration is small, the vicinity of the percola- 
tion fixed point induces a decrease of 1 jv below its 
expected disordered value. Indeed, the percolation 
fixed point is characterized by \ jv ~ 1.12 [7]. 



4. Conclusion 

We have presented, from a numerical study, the 
influence of bond dilution on the critical properties 
of the 3D Ising model. The universality class of 
the disordered model is modified by disorder but 
its precise characterization is difficult because of 
the competition between the different fixed points 
which induce crossover effects, even for relatively 
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Fig. 5. Effective exponents (l/^) e £f as a function of 1/L m j n 
for p = 0.95, 0.9, 0.8, 0.75, 0.65, 0.6, 0.45 and 0.36. The 
error bars correspond to the standard deviations of the 
power-law fits. The arrows indicate the values of 1/v for 
the pure model [5] and the site-diluted one [3], 

large lattice sizes. The next step is to accurately 
locate the best dilution which minimizes these 
crossover effects. 
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